Structure of Growing Networks: Exact Solution of the Barabasi— Albert's Model 
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We generalize the Barabasi-Albert's model of growing networks accounting for initial properties 
of sites and find exactly the distribution of connectivities of the network P{q) and the averaged 
connectivity q{s,t) of a site s in the instant t (one site is added per unit of time). At long times 
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at g ^ oo and q{s,t) ~ (s/t) ^ at s/t — > 0, where the exponent 7 varies from 2 to 00 



depending on the initial attractiveness of sites. We show that the relation /3(7 - 
exponents is universal. 
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Recently opened intriguing scaling properties of a vari- 
ety of growing networks (WWW, neural, social, citations 
of scientific papers, etc., see [p]-p^ and references therein) 
challenge us to find general reasons of such behavior. In 
fact, the studying of networks has a long story |14-1^] 
but only now we become feel ourselves as subjects inside 
of the world of developing networks and start appreciate 
significance of the evolving network problem ^ 6| |2C[] . 

Mostly, the interest is concentrated at the distribution 
of shortest pathes between the different sites of a network 
101 and at the distribution of the number of connections 
with a site P{k) [||-^. The second quantity is obviously 
simpler to obtain but even for it, in the case of the net- 
works with scaling behavior, no exact results are known. 

A clear model of a growing networks was proposed by 
Barabasi and Albert The model presents the most 

simple mechanism of self-organization of a growing net- 
work into a free-scale structure. At each time step a new 
site is added and is connected with old sites of the net- 
work through m new links. The probability of an old 
site to get a new link is proportional to the total num- 
ber of connections with this site. It was found in 
that P{k) ~ k~'' at long times and large connectivities, 
where 7 = 8 (the approach that seemed at first sight be 
approximate - "mean field" - was used). This value is 
close to that one observed in the network of citations , 
but other examples show different values of 7. Introduc- 
tion of the aging of sites changes 7 and may even 
break the scaling behavior [|l2|,^ . 

In the present letter, we generalize the Barabasi- 
Albert's model in natural way and solve it exactly, i. e. 
we find the exact form of the distribution of connectivi- 
ties and some other related quantities of the network. We 
show also that the appearing scaling relations are valid 
not only for this particular model but for a wide class of 
growing networks. 

Let us formulate the model. At each time step the 
new site with m connections with some old sites appears. 
Hence, in each instant t = 1,2,..., the network consists 
of t sites connected by mt directed links. Note that we 
consider not a total number of connections with a 
site but Qs, the number of incoming links, where s is the 



number of a site. Therefore, for brevity, we call it the 
connectivity of a site but not kg as usually. One intro- 
duces the following rule to distribute m new links among 
t old sites in the instant t + 1. Let the probability that 
a new link points to a given old site s be proportional to 
the following characteristic of the site: 
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that may be named its attractiveness. All sites are borned 
with some initial attractiveness A^^') but afterwards it in- 
creases because of the qs term. In the particular case of 
the Barabasi-Albert's model, ^'"^ = m. 

In fact we consider the following general problem. In 
each instant, m new particles (i. e. incoming links) have 
be distributed between an increasing number (by one per 
time step) of boxes (i. e. sites) acording to the intro- 
duced rule. One has to emphasize, that in our formula- 
tion, there is no difference, from which particular site - 
old or new or even from outside of the network - these 
new links come. Indeed, we are interested only in the 
statistics of incoming links. Therefore, new links may 
also appear between old sides. 

It is convenient to assume that in the initial instant t — 
1 we have only one site with the connectivity m. Then 
the total attractiveness of the network in the instant t 
is A^ = {m + A^°'>)t = (1 -I- a)m, where a = A'^°'> /m. 
The resulting behavior at long times is independent of 
the initial condition. 

Let us derive the equation for the distribution p{q, s, t) 
of connectivities q at the site s in the instant t. The prob- 
ability that a new link is connected with the site s is Wg = 
As /As. The probability for the site to receive exactly I 

new links of m injected is = (7)PFi (1 - Wg)"'~\ 

Note that we allow multiple links unlike the original 
Barabasi-Albert's model, i. e. the connectivity of a given 
site may increase by more than one in the same instant. 
That is inessential at long times, when probability to re- 
ceive more than one link of m introduced simultaneously 
is vanishingly low. Hence, the connectivity distribution 
of a site obeys the following master equation: 
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Eq. is supplied with the initial condition p (q, s, s) = S (q), which means that sites are horned with zero connectivity 
(i. e. without incoming links in our definition). 

The distribution function of connectivities in the all network is 



1 

P(g,t) = -^p(g,w,0 . 

u=l 

Summing up Eq. over s from 1 to t, one gets 

{t+l)P{q,t + \)-p{q,t+\,t+l)^-[t-\^]P{q,t) + ^— P(g-l,i) + - 

\ 1 + a J 1 + a \t 

In the long-time limit, after the transition to continuous-time approximation, we obtain 
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(1 + a) t— (q, t) + {l + a)P {q, t) + (q + am) P {q,t) - {q - I + am) P {q - l,t) ^ (l + a) S (q) . (5) 
ot 

Finally, assuming that the limit, P (q) = P{q,t^oo), exists, we get the following equation for the stationary 
connectivity distribution: 



(1 + a) P (q) + iq + ma) P (q) - {q - 1 + ma) P (q - 1) = {1 + a) 6 (q) 



(6) 



To solve Eq.(§) one may use Z-transform of the distri- 
bution function: 



Then one gets from Eq. 



$(z) = ^P(g), 

9=0 



(7) 



following form: 



z{l - z) — +ma{l - z)^ + {1 + a)^ ^ 1 + a. (8) 
dz 

The solution of Eq. (g) that is analytic at z = has the 
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2Fi[l,ma;2+ {m + l)a;z], (9) 



where 2-F1 [ ] is the hypergeometric function. Using its 
expansion ||22| in z, we obtain, comparing with Eq.(|^), 



P(g) = (l + a) 



r [(m -I- 1) a + 1] T{q + ma) 

r(ma) T[q + 2 + {m + l)a] 



(10) 



that is our main result (see Fig. 1). In particular, when 
a = 1 , that corresponds to the case As = m + q,, = kg , 
studied in we get exactly 



p{q) = 



2m (to -t- 1) 



iq + m){q + m + l){q + m + 2) 



(11) 



This expression in the limit q 00 approaches the cor- 
responding result of obtained in the frames of an 
approximate scheme, but the prefactors are different. 
In fact, the "mean field" approximation, used in jsj, is 
equivalent to continuous-g approximation in our discrete- 
difference equations. Indeed, if we replace the finite dif- 
ference with a derivative over q, we get the expression 
obtained in 0J|]. 



At TOO + q ^ 1, the distribution function ( p^ takes 
the form: 
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P(g)-(l + a) 



r[(m+l)a+l] 



(<7 + ma)-('+"^ . (12) 



where A*^") is the initial attractiveness of a site. 



T (ma) \i / \ ' j^^^ gj^j distribution function p(q,s,t) of con- 

Therefore, we find the scahng exponent 7 of the distri- nectivities at the site s. At long time_s, t » 1, keeping 



bution function: 

7 = 2 + a 2 + ^("V^ 



only two leading terms in 1/t in Eq. (y), one gets 
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p{q,s,t+ 1) 
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P + il + a)t P (9 - 1, ^) + O j 
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(1 + a)i 

Assuming the scale of time variation to be much larger then 1, wc can replace the finite ^-difference with a derivative 



dp 

(1 + a) t— {q, s, t) = (g - 1 + am) p {q - 1, s,t) - {q + am) p {q, s, t) 



(15) 



Finally, using Z-transform in the similar way as before, 
we obtain the solution of Eq. (p^), i. e. the connectivity 
distribution of individual sites: 
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Hence, this distribution has an exponentional tail. Now 
one may get also the expression for the average connec- 
tivity of a given site: 
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q{s,t) = ^qp{q,s,t) 
9=0 



-l/(l+a) 



- 1 



(17) 



Thus, at a fixed time t the average connectivity of an 



old site s <^ t depends upon its age as 



^ , where 



the exponent (3 = \/ {1 + a). Therefore, we have the fol- 
lowing relation between the exponents of the considered 
network: 



/3(7-l)-l 



(18) 



that was obtained in |2l[ with some particular form of 
p{q,s,t). 

We can show that Eq. (^Sj) is universal and may be ob- 
tained from the most general suggestions. In fact, we as- 
sume only that the averaged connectivity q{s,t) and the 
distribution of connectivities P{q) show scaling behavior. 
Then, in the scaling region, the quantity of interest, i. e. 
the probability p{q, s, i), has to be of the following form: 
p{q,s,t) = {s/t)^^f{q{s/t)^^). Here, Ai = A2 because 
of the normalization condition for p{q, s, t) at a fixed s, 
X]^o P ~ ^- Then, the relation g(s, t) oc {s/t)~^ 

leads to Ai = A2 — [3 (we use the definition (|l7|)), and 
finally, inserting p{q, s, t) in such a form into the relation 
P{q) oc at large q and t, one gets the relation ([l^). 

The particular form of the scaling function /{£,),£, = 
q{s/t)~^ depends on the specific model of the growing 
network. In the case under consideration, it follows from 
Eq. (H) that 
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at s/t 0, q ^ 00 and the fixed q{s/t)f^ . (We used here 
that T{am + q)lq\ q"™-! at g ^ 00.) 

In the limit of zero initial attractivity of sites (a = 0) 
all new sites connect only with the first one, since all 
other sites have no chance to get a new link. In this case, 
Eqs. (|l|) and (|l|) give 7 = 2 and /3 = 1. For a = 1, i. 
e. for the Barabasi-Albert's model, 7 = 3 and (3=1/2 
Finally, when a ^ 00, i. e. all sites have equal at- 
tractivity all time, and the scaling breaks, one sees from 
Eqs. (|l|) and (^ that 7 ^ 00 and /? ^ 0. Note that 
the ranges of variation of 7, (2, 00), and (3, (1, 0), are the 
same as for the network with aging . 

We see that the "approximate" approach of gives 
the proper values for the critical exponents (see also |2|] 
for the network with aging). Is that only an occasional 
coincidence? In fact, in these papers, some particular 
form of p{q,s,t) was used (i) to derive the equation for 
the averaged connectivity and, therefore, for (3, and (ii) 
to obtain relation between 7 and (3. One can noticed, 
however, that (i) the equation for the averaged connec- 
tivity may be obtained for arbitrary p(g, s, t) and (ii) the 
relation between the scaling exponents is universal, as we 
have shown above. Thus, the analytical results for the 
scaling exponents obtained in ||^,||,^ turn to be exact. 

A two-parameter fitting was proposed in |Q to de- 
scribe the observed distribution of citations of scientific 
papers ||^. One sees, that the connectivity distribution 
of the considered growing networks is of the quite differ- 
ent form. It seems, that the difference occurs because 
we study the growing structure unlike the approach |2|] 
although the question is open. 

In conclusion, we have found the exact form of the 
connectivity distribution and the averaged connectivity 
of sites of the natural generalization for the Barabasi- 
Albert's model. The considered growing network is self- 
organized into the free-scale structure. The input ffow of 
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the new links is distributed between the increasing num- 
ber of sites. In fact we have considered the new version 
of a sandpile problem. The scaling exponents are deter- 
mined by the value of the initial attractiveness ascribed 
to every new site. Depending on this quantity, the scaling 
exponent 7 of the connectivity distribution takes values 
from 2 to cx). We have shown that the relation between 
the scaling exponents 7 and (3 is valid for a wide class 
of growing networks. Our simple solution of this prob- 
lem lets hope to obtain other exact results for growing 
networks. 

The following question remains open. We have found 
that the exponent 7 varies between 2 and 00. Why is 7 
only in the range of 2-3 in real networks fip^? 
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FIG. 1. Log-log plot of the distribution of incoming links 
of sites, m = 1. 1) a = 0.001, 2) a = 0.05, 3) a = 1.0 (the 
Barabasi-Albert's model), 4) a = 2.0, 5) a = 4.0. 
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